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Abstract 


A 6-ha field at Aberystwyth, UK, was converted in 2012 from semi-improved grassland to Miscanthus x giganteus 
for biomass production; results from transition to the end of the first 3 years are presented here. An eddy covari- 
ance sensor mast was established from year one with a second mast added from year two, improving coverage 
and providing replicated measurements of CO, exchange between the ecosystem and atmosphere. Using a sim- 
ple mass balance approach, above-ground and below-ground biomass production are combined with partitioned 
CO, fluxes to estimate short-term carbon deltas across individual years. Years one and two both ended with the 
site as a net source of carbon following cultivation disturbances, cumulative NEE by the end of year two was 
138.57 + 16.91 g Cm ~. The site became a cumulative net sink for carbon by the end of June in the third 
growing season and remained so for the rest of that year; NEE by the end of year three was —616.52 + 
39.39 g C m *. Carbon gains were primarily found in biomass pools, and SOC losses were limited to years one 
(-1.43 Mg C ha! yr) and two (-3.75 Mg Cha! yr’). Year three saw recoupment of soil carbon at 
0.74 Mg C ha“! yr! with a further estimate of 0.78 Mg C ha“! incorporated through litter inputs over the 
3 years, suggesting a net loss of SOC at 3.7 Mg ha“! from a 0- to 30-cm baseline of 78.61 + 3.28 Mg ha~', down 
4.7%. Assuming this sequestration rate as a minimum would suggest replacement of cultivation losses of SOC 
by year 8 of a potential 15- to 20-year crop. Potential coal replacement per hectare of harvest over the three-year 
study would offset 6-8 Mg of carbon emission, more than double the SOC losses. 
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Introduction 


There is much debate surrounding the potential for car- 
bon sequestration into soils under perennial energy 
crops such as Miscanthus x giganteus (hereafter Miscant- 
hus), particularly when planted into agricultural grass- 
lands (McCalmont et al., 2015). Some studies find 
increased soil organic carbon (SOC) when comparing 
Miscanthus plantations to grassland (Hansen et al., 2004; 
Clifton-Brown et al., 2007; Schneckenberger & Kuzya- 
kov, 2007); others find SOC to be unchanged or even 
reduced (Zimmermann et al., 2012; Zatta et al., 2013). 
Regardless of direction, the significance of any changes, 
particularly short term, is often difficult to demonstrate 
due to small changes in large volumes and the limita- 
tions of direct sampling (Smith, 2004; Kravchenko & 
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Robertson, 2011). These challenges are typically com- 
pounded by a reliance on adjacent land taken to repre- 
sent baseline conditions and incompatibility of 
sampling techniques and depths in comparisons. Sam- 
pling to a fixed depth when comparing Miscanthus to 
grassland or following cultivation can be criticized for 
not accommodating changes in soil bulk density which 
can exaggerate changes in carbon stocks (Gifford & 
Roderick, 2003; Rowe et al., 2015), although techniques 
such as equivalent soil mass (ESM) are increasingly 
employed (Ferchaud et al., 2015; Richter et al., 2015; 
Rowe et al., 2015). Understanding the source of soil car- 
bon can be more reliable when investigating a change 
from a conventional C3 to a C4 crop, such as Miscant- 
hus, as the C4 vegetation shows less discrimination 
against the naturally occurring °C isotope in 
atmospheric CO, (Balesdent et al., 1987). Soil carbon 
derived from Miscanthus shows less depletion of this 
isotope when compared to an atmospheric standard in 
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comparison with C3-derived carbon. Studies make it 
apparent that there is an initial and fairly rapid decom- 
position of soil carbon following cultivation disturbance 
and soil respiration priming (Cheng et al., 2003; Kuzya- 
kov, 2010; Hopkins et al., 2013) while isotopic analyses 
show that this can be recouped to varying extents over 
the coming years by C4 inputs from Miscanthus (Hansen 
et al., 2004; Dondini et al., 2009; Zatta et al., 2013). Com- 
plete replacement or net sequestration will occur only 
where C4 inputs outpace C3 losses and, despite being 
generally reported as a mean rate over several years, 
recoupment rates will not be constant. Agostini et al. 
(2015) stress that soil carbon input rates are a function 
of crop yield and time, increasing as litter drop and root 
turnover increase with maturing crops and accumulate 
in litter and biomass pools. Direct sampling from infre- 
quent, limited soil coring is less capable of capturing 
these short-term dynamics, so there is a need for a more 
sensitive approach to estimate dynamic changes in soil 
carbon. Results might also indicate a minimum crop 
lifetime where carbon gains in biomass might outweigh 
any recycling of soil carbon to the atmosphere or vege- 
tation. Previous studies in forest systems (Raich & 
Nadelhoffer, 1989; Giardina & Ryan, 2002) have sug- 
gested a mass balance approach to estimating below- 
ground carbon allocation where soil respiration is 
related to estimates of total carbon input to the soil 
through litter drop, harvest residue and root turnover. 
Smith et al. (2010) propose taking advantage of eddy 
covariance measurements of whole-system carbon 
exchange (Baldocchi et al., 1988) and comprehensive 
assessments of biomass production in croplands to pro- 
duce a net ecosystem carbon budget. Smith et al. (2010) 
suggest that total photosynthetic uptake (GPP) should 
equal ecosystem respiration and total biomass (above 
ground and below ground) production (NPP) in a stable 
system, with some allowance for losses through soil 
water leaching or erosion. However, agricultural sys- 
tems are rarely stable over time in terms of soil carbon 
or organic matter, particularly following cultivations. 
Despite ploughing in of previous crop residues, con- 
stant stimulation and release of soil carbon through dis- 
turbance primed respiration can result in net carbon 
release from soils to atmosphere and soils under annual 
arable cropping typically show lower carbon stocks 
when compared to perennial systems such as grasslands 
or Miscanthus (Felten & Emmerling, 2012; Poeplau & 
Don, 2014). 

In this present study, we assume that discrepancies 
between the amount of carbon taken in by the vegeta- 
tion and biomass production plus ecosystem respiration 
can be taken to indicate losses or gains to soil carbon, 
that is if it appears biomass production and respiration 
together outweigh photosynthetic uptake from the 


atmosphere above the canopy then there must be a net 
loss of carbon from the soil. While immaterial to this 
mass balance approach, it should also be considered 
that the growing crop may directly recycle this respired 
soil carbon into plant biomass, that is vegetation utiliz- 
ing below canopy availability of CO) respired through 
decomposition of soil carbon. Buchmann & Ehleringer 
(1998) demonstrated this uptake of soil respired CO, 
directly into C3 and C4 crop leaves. They showed pho- 
tosynthetic demand occasionally outstripping total soil 
respiration and depleting subcanopy CO, levels below 
tropospheric concentrations while Brooks et al. (1997) 
reported greater subcanopy diffusion gradients and 
understory leaves deriving 5-6% of their carbon from 
respired CO}. 

In addition to a mass balance estimation of annual 
carbon stocks, eddy covariance allows real-time assess- 
ment of source/sink dynamics and interyear compar- 
isons of the impact of a developing perennial crop on 
atmosphere/ecosystem carbon exchange. In this study, 
we present the results of three-year measurements using 
eddy covariance and biomass sampling at a commer- 
cial-scale land-use change experiment in Wales, UK, 
from semi-improved agricultural grassland to Miscant- 
hus. The study site is one of five flagship land-use con- 
version sites developed across the UK within the wider 
Carbo-Biocrop (http://www.carbo-biocrop.ac.uk) and 
ELUM (http://www.elum.ac.uk) projects with the pri- 
mary aim of providing validation data from the first 
2 years to improve greenhouse gas flux estimation in 
the ECOSSE model (Bell et al., 2012; Harris et al., 2014). 
In addition to eddy covariance and carbon stock assess- 
ments, trace gas fluxes (N?O and CH4) were sampled 
monthly across all sites from closed static chambers 
across all five sites, the results of model validation 
against these data were presented in Dondini et al. 
(2015). Despite the completion of these projects at the 
end of 2013 after 2 years of measurements, the site con- 
tinues to be monitored for carbon fluxes with results 
from the first 3 years presented here. 


Materials and methods 


Site description 


The study site is a 7.41-ha field at Penglais, Aberystwyth, in 
mid-Wales (52°25/17" N 4°04'14" W); ~ 110 m a.s.l. with the 
western edge of the site sitting at the top of a steep decline 
down to the coast. Soil type is a sandy loam, dystric cambisol 
over Denbigh series bedrock with a mean pH of 5.9. Climate is 
temperate with 30-year local annual averages of 158 days with 
rain, 1074.7 mm total rainfall and max/min temperatures of 
13.5/6.7 °C [Gogerddan 1981-2010 averages (www.metoffice. 
gov.uk)]. Prior land use was semi-improved grassland which 
was ploughed and resown regularly over the previous 
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30 years. The last time was 6 years prior to this experiment, 
which meant the site was due for recultivation at this point 
with only the novel crop type straying from what would have 
been normal practice anyway. Due to constraints of homogene- 
ity imposed by the eddy covariance (EC) technique, the land- 
use change (LUC) area was restricted to a central 5.71 ha with 
two control areas of perennial ryegrass-dominated grassland 
(RGE & RGW) totalling 1.59 ha retained at either end for com- 
parisons of biomass production (Fig. 1). For baseline, the LUC 
area was considered in two blocks (LUCE & LUCW) as this 
related directly to the measurement footprints of the individual 
EC masts (Fig. 5) with differences in soil parameters and bio- 
mass production between Miscanthus and grassland areas 
investigated using anova with block differences separated 
using Tukey HSD analysis. 


Land-use conversion 


The existing grass was intensively grazed across the entire site 
to ground level by sheep at 10 head ha~' by the 20 February 
2012; the sheep were then removed and the land-use change 
area (LUC) sprayed with glyphosate at 1.5 kg ha™' on the 16 
March 2012, killing off the remaining sward. This area was 
ploughed to ~ 20 cm depth on the 4 April 2012, incorporating 
the dead grass and roots into the soil, followed by power har- 
rowing on the 23rd April. Miscanthus x giganteus rhizomes were 
commercially planted on the 24th April by International Energy 
Crops (Market Drayton, UK) at a target density of 
~ 16 000 ha~'. Table 1 shows a list of herbicide application in 
the LUC area over the 3 years, no fertilizers were applied to 
either land use. 


Soil texture 


Baseline soil parameters were determined by core sampling 
prior to cultivation. Twenty cores (1=5 per block) were 
extracted across the site using a noncompressive 
1m x 0.085 m diameter, tractor-driven soil column cylinder 
auger (Eijkelkamp, Giesbeek, The Netherlands). Cores were 
separated into 0.15-m sections, dried at 60 °C until constant 
weight, crushed and sieved to <2 mm to remove stones. Bulk 
densities with and without stone were calculated for each 15- 
cm section. Finer texture analysis was carried out by laser 


Table 1 Herbicide application across the three-year study 
period. No further pesticides or fertilizer were applied in any year 


Brand Application 
Date Chemical name rate (ha!) 


16-03-2012 Glyphosate Round-up 1.5 kg 

17-05-2012 Isoxaflutole 100 g kg" Cadou 0.75 kg 
Flufenacet 480 g kg! Star 

26-07-2012 Metsulfuron-methyl Jubilee SX 30 g 
Bromoxynil/ioxynil Oxytril IL 
Fluroxypyr Starane I 0.5 L 
as 1-methyl heptyl 

03-04-2013 Glyphosate Round-up 2 kg 
Chlortoluron 620 g lt™™ Steel 3L 
Diflufenican 22.5 g lt~! 

16-05-2014 Pyroxsulam Broadway 0.25 kg 
70.8 g kg! Star 
Florasulam 14.2 g kg~! 


ELLE 


150 200 250 Meters 


Fig. 1 Experimental layout of the Aberystwyth site, LUCW/E indicates the land use area converted to Miscanthus in 2012 
(W = west, E = east) while RGW/E indicates retained original grassland. Open squares show biomass sampling quadrats for 2012. 
EC1 and EC2 show eddy covariance mast locations with the measurable fetch boundary line demarcating the ideal flux footprint geo- 
referenced for EC data quality control. The meteorological station can be seen near the centre of the field and power supply tower at 


the eastern edge. 
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diffraction particle size analysis (Beckman and Coulter LS200, 
Beckman Coulter Ltd, High Wycombe, UK) on a total of 24 
subsamples: 12 locations at two depths [0-0.15 m and 0.15- 
0.30 m (n = 3 at 2 depths per block)]. Subsamples from the two 
depths (n = 3 at 2 depths per block) were analysed for soil 
organic carbon (SOC) through oxidation with potassium 
dichromate to establish a baseline carbon content to 0.3 m. 


Baseline biomass 


Intensive grazing reduced the grass crop to ground level at the 
site before conversion leaving minimal above-ground biomass. 
To assess total baseline biomass, square sections were cut 
through into the top 0.15 m of soil which captured this and the 
below-ground biomass together. Ten of these were taken with 
dimensions of 0.10 x 0.10 m x 0.15 m deep, randomly dis- 
tributed across the study site. Samples were then washed clean 
of soil and dried to constant weight. Dry matter contents for all 
biomass collections were assessed by oven drying at 60 °C until 
constant weight. 


Miscanthus yield and below-ground biomass (LUC) 


For years one and two, eight randomly distributed quadrats 
(n = 4 per block) were established within the LUC area at the 
beginning of each growing season (Fig. 1), for these first 
2 years significant differences in biomass production were 
assessed between the Miscanthus blocks (LUCW/E) and 
between these and grass production on the perennial rye grass 
(RGW/E) blocks. For year three, due to limited resources fol- 
lowing the end of the ELUM project, quadrat numbers were 
reduced to five within the LUC area and differences between 
blocks or retained grassland were no longer assessed. Quadrats 
were 2 m long and 1.22 m wide (2.44 m°) covering two rows of 
Miscanthus and two inter-rows, and on average, there were six 
plants captured in each quadrat. These plants were repeat mea- 
sured for canopy height at approximately weekly intervals (see 
Fig. 4) and destructively harvested following senescence at the 
end of each growing season to determine peak yield. The year 
one crop was cut and left in the field, so all production was 
considered to be litter drop into the following year. A second 
assessment of biomass from eight further quadrats was carried 
out in early March of the following year to determine spring 
harvestable yield and overwinter litter drop. The crop from the 
year two growing season was commercially harvested with no 
direct measurement of offtake beyond this subsampling. The 
year three crop was also commercially harvested in early 
spring of the following year, but all material removed from the 
site was this time weighed directly. 

Differences between peak and harvest yield are taken to rep- 
resent overwinter litter drop, primarily leaf loss. At the end of 
year three, litter stock remaining on site was assessed from 20 
randomly located quadrats (0.26 m~’), differences between this 
and the estimate of total litter drop over the 3 years were taken 
to represent the quantity that had decomposed during the 
three-year period. 

Rhizome and root development were determined from each 
of the sampling quadrats during the dormant period between 


growing seasons. A section of soil, one from each quadrat, was 
removed centred on one Miscanthus plant with two associated 
inter-rows. Plants were randomly selected from within each 
quadrat following removal of the above-ground biomass. This 
produced eight root cores 0.5 m° and 0.30 m deep, capturing 
entire rhizomes and the visible coarse roots. Baseline rhizome 
biomass was determined by dividing the known fresh bulk 
weight of the planted material by the planted area with an 
assessment of moisture content carried out on ten 100 g sub- 
samples. 


Scaling from plant to field scale 


Scaling from individual plants to results per hectare can be 
problematic in a commercial-scale row crop due to gaps within 
and between rows. To provide a scaling parameter of plants 
ha™!, five 20 x 20m quadrats were marked out at random 
locations across the LUC area and all Miscanthus plants con- 
tained were counted to provide an estimate of plants ha~' from 
400 m?. This was carried out for the first 2 years to accommo- 
date any overwinter losses during establishment. 


Retained grassland yield and below-ground biomass 
(RG) 


Biomass production on the control areas of grassland (RGW/E) 
was monitored for the first 2 years of the study. Steel hoops 
with an area of 0.26 m? were used as randomly distributed 
sample quadrats; all grass contained within them was cut to 
ground level and removed for drying; five of these samples 
were taken from each of the retained grassland areas. This was 
carried out twice during each year, once in June immediately 
before the grass was mown and removed (no grazing was car- 
ried out during the study period), and again at the end of the 
growing season. The sum of these weights represents cumula- 
tive above-ground biomass production of these retained grass- 
land areas from a baseline of essentially zero when the sheep 
were finally removed. Below-ground root biomass in the RG 
areas was sampled at the same time as the Miscanthus from 
four randomly distributed cores taken within each block, 
0.16 m* x 0.15 m deep. 


Instrumentation 


Meteorology. A meteorological station was established at the 
centre of the fetch area from January 2012 with data recorded 
to a CR1000 datalogger (Campbell Scientific Inc. (CSI), Logan, 
UT, USA). Measurements included net radiation (Rn, W m7”), 
NR Lite net radiometer (Kipp and Zonen, Delft, Netherlands); 
photosynthetic photon flux density (PPFD, umol ms), 
SK215 quantum sensor (Skye systems, Llandrindod Wells, UK); 
precipitation (P, mm), Young’s 52203 tipping bucket rain gauge 
(R.M. Young, Michigan, USA); soil temperature (Ta, °C), TCAV 
(CSI); and volumetric water content (vwc, %) CS616 sensors 
(CSI). Soil moisture and temperature were measured at three 
locations orthogonal to each other five metres out from the sen- 
sor mast into the field. Two of these measured moisture and 
temperature at two depths, 0.025 and 0.25 m, the deeper probes 
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being immediately above a layer of gravel in the shallow soil. 
Also located centrally were soil energy balance sensors, these 
consisted of a horizontal CS616 soil moisture probe with its 
upper tang at 0.025 m below the soil surface with TCAV tem- 
perature probes 0.5 m either side and HFP01SC (CSI) heat flux 
plates below these at 0.08 m. 


Eddy covariance 


Fetch. The area of the site that satisfied eddy covariance 
topographical assumptions of minimal slope and maximum 
downwind fetch availability restricted the ideal EC fetch to 
the flattest 3.9 ha of the LUC area (Fig. 1). The field itself 
was chosen for its alignment lengthwise to the prevailing 
south westerly winds; to maximize this the mast was 
installed at the eastern edge of the field (EC1 in Fig. 1). 
Lengthwise fetch in this prevailing wind direction was 255 m 
with 80 m to either side. In practice, the output of the foot- 
print model showed that measurements from the south-west 
rarely extended beyond the LUCE block boundary (Fig. 5) 
leaving LUCW largely unrepresented in 2012. A second iden- 
tical eddy covariance mast was installed in January 2013 at 
the opposite end of the fetch (EC2 in Fig. 1), distance 
between the two masts was 180 m. The second tower signifi- 
cantly improved spatial and temporal coverage over the 
maturing crop and took advantage of a regular nocturnal off- 
shore reversal of wind direction and accompanying good 
night-time turbulence. This twin tower set-up also offered a 
rare opportunity in eddy covariance studies to directly cap- 
ture sampling uncertainty. 


EC sensors 


The eddy covariance systems (EC150/CSAT3A OPEC system, 
CSI) consisted of a CSAT-3A sonic anemometer and EC150 
infrared gas analyser with air temperature (Ta °C) and rela- 
tive humidity (RH, %) monitored by an HMP155A (CSI) 
along with soil surface moisture and temperature five metres 
forward of each flux mast at a depth of 0.025 m using CS616 
and TCAV sensors (CSI); data were recorded to a CR3000 
datalogger at 20 Hz and later processed to 30-min averaging 
intervals. Sensor heights needed adjustment to accommodate 
the growing crop, allowing them to remain high enough 
above the canopy to avoid the roughness sublayer of turbu- 
lence where fluxes would be difficult to resolve but low 
enough to restrict the fetch to within the crop boundaries. A 
target above canopy sensor height was therefore set at 2 m, 
to allow adjustment the EC sensors were mounted on a mov- 
able cross arm mounted onto the central 6-m spar. Crop mea- 
surements were taken approximately weekly during the 
growing season and when canopy heights indicated this dis- 
tance had narrowed to below 2-m sensors were raised 
accordingly. CO, storage in the canopy during periods of low 
atmospheric turbulence was considered negligible over time 
and only above-canopy NEE was measured, with the stan- 
dard sign convention that negative values represent uptake of 
CO, from the atmosphere and vice versa. Mean flux rates 
of CO» (umol CO2 ms!) are converted to mass sums of 


carbon (g CO-C m~? hh™') for budget integrations. The 
EC150s were calibrated as new from the factory on installa- 
tion and subsequently zero and span calibrated at the begin- 
ning and end of each growing season using zero-grade air 
and a 500 ppm CO, standard (BOC speciality gases, Surrey, 
UK), drift was minimal with an average gain across both 
towers of 0.003 ppm. Sensors were reset at each calibration, 
and results are not corrected for the minimal gain. During 
the winter between the second and third years, the gas analy- 
sers were returned one at a time for factory calibration. 
Power was supplied by wind (Rutland 913 wind charger, 
Marlec Engineering, Northants, UK) and solar power (80W 
solar PV, RS Components Ltd, Northants, UK) mounted on a 
5-m scaffold tower 20 m downwind of EC1 (Fig. 1). 


EC data quality control 


The 20 Hz data were processed into 30-min average flux rates 
(Baldocchi, 2003) using EDDYPRO software (EddyPro® version 
4.2.0, LI-COR bioscience, Lincoln, NE, USA). Quality control 
flagging policy (0, 1, and 2) for each half-hour average fol- 
lowed Mauder & Foken (2004); spike detection followed Vick- 
ers & Mahrt (1997); vertical wind speed (W) retained at <5 
standard deviations (a) from the mean over the half hour, all 
other variables at 3c. Detrending of turbulence fluctuations 
used block averaging. Crop canopy dynamics and associated 
sensor heights are included in the processing by automated 
reading of metadata files. Coordinate rotation was by double 
rotation rather than planar fit (Wilczak et al., 2001), which 
would require undisturbed sensor alignments over longer 
timescales to produce reference wind vector parameters. Den- 
sity fluctuations were corrected with the WPL term (Webb 
et al., 1980; Leuning, 2007). Co-spectral analysis and correc- 
tion of low- and high-pass filtering effects were carried out 
following Moncrieff et al. (1997, 2004). Cross-wind corrections 
were handled internally by the CSAT-3A, velocity bias in the 
anemometer calibration was assessed before installation with 
deviations from zero included in the processing calculations. 
This half-hourly data set was then further quality controlled 
for site-specific parameters with values flagged at 2 rejected 
immediately. Sensor obstruction, for example precipitation, is 
indicated by mean signal strength between transducers, 
return signals <85% of the outgoing were rejected. Flux foot- 
print was estimated using the models of Kljun et al. (2004) or 
Kormann & Meixner (2001) where turbulent friction velocity 
was <0.2 m s™!. cis software (ArcMap 10.1 ESRI, Redlands, 
CA, USA) was used to determine distance to the outer limit 
of the acceptable fetch within 36 ten-degree increments; com- 
pletely unacceptable directions, due to flow distortion by the 
sensor masts themselves or interference from surrounding 
trees/hedgerows, were excluded immediately; remaining 
increments were acceptable where 70% of the measured flux 
was determined by the footprint model to have come from 
within the acceptable fetch for each ten-degree increment 
(Fig. 5). Post processing data handling and quality control 
were carried out using script programming in the R statistical 
language (R Foundation for Statistical Computing, Vienna, 
Austria). 
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Gapfilling 


Data set gaps created either through original data losses or 
rejected through quality control need to be filled for estima- 
tions of long-term carbon budgets. As can be seen in Fig. 5, 
where wind directions excluded data collected at EC1 they 
were often acceptable at EC2, in these instances and others 
where quality control has rejected data from EC1 but not EC2 
this data set was used to fill gaps in EC1. The remaining gap- 
filling, along with flux partitioning of NEE into ecosystem res- 
piration (Reco) and gross primary production (GPP), was 
carried out using the FLUXNET standard online gapfilling tool: 
http://www.bgc-jena.mpg.de/bgi/index.php/Services /REd- 
dyProcWeb. 

Underestimation of fluxes during periods of low turbulence 
was avoided by filtering data under derived friction velocity 
(u*) thresholds (Reichstein et al., 2005). Energy balance closure 
was considered in ordinary least square regression between the 
latent plus sensible heat (LE+H) and net radiation minus soil 
heat storage (Rn—G) using all half-hour values retained in each 
year (Wilson et al., 2002). Canopy storage of heat or carbon 
dioxide was not measured and is not included. 


EC data set comparison 


The footprints of each sensor mast were predominantly inde- 
pendent of each other, particularly during westerly wind direc- 
tions, offering direct comparisons between masts so estimating 
sampling uncertainty across the site [see Hollinger et al. 
(2004)]. Site- and mast-specific quality control was applied to 
both the EC1 and EC2 data sets, including filtering values 
recorded where u* was below the site-specific threshold, before 
appropriate paired measurements were selected out from the 
retained data. These were measurements acceptable from both 
towers at the same half-hour time points and thereby driven by 
near identical climatic factors, these resulting data pairs were 
compared using linear regression. 


Carbon budget 


For a spring harvested perennial crop such as Miscanthus, a 
year running from January to December may not be the most 
appropriate time span to consider interannual comparisons of 
net carbon exchange or budget sums. For these results, periods 
running from the beginning of March in 1 year to the end of 
February in the next (nominally harvest to harvest) are consid- 
ered. Therefore, we present results from year one (17/3/2012 
to 29/2/2013), year two (1/3/2013 to 28/2/2014) and year 
three (1/3/2014 to 28/2/2015). Year one begins 17th March as 
this is the day after spraying off the original grassland at the 
beginning of the land-use change. 

A mass balance approach was taken to estimate changes in 
below-ground carbon stocks, both for total below-ground car- 
bon including biomass (TOC) and for soil organic carbon 
(SOC). Conceptually, if respiration and the carbon content of 
biomass production (above ground and below ground) are sub- 
tracted from total photosynthetic uptake then SOC losses or 
gains can be estimated. Negative values would indicate a loss 


of soil carbon, positive would indicate a gain (note this is now 
opposite to the eddy covariance convention). This simple con- 
cept yields two equations. .. 


ATOC = (1) 


= -— (Ra z ANPP)) + *) : 


(ATOC — (BNPP + rz)) 
t 


ASOC = (2) 


where 

ATOC = change in total below-ground organic carbon (in- 
cluding biomass) 

GPP = gross primary productivity (total CO, uptake in pho- 
tosynthesis) 

Reco = total ecosystem respiration 

ASOC = change in soil organic carbon 

a = any vegetation added to the below-ground carbon pool 
in cultivation, that is original grass killed and ploughed in, 
subsequent weed roots killed with herbicide control and addi- 
tions through rhizome planting and litter decomposition 

ANPP = above-ground net primary productivity (biomass 
production) 

BNPP = below-ground net primary productivity (root/rhi- 
zome biomass gain) 

rz = planted rhizome 

t = time span 

é = losses through soil erosion and dissolved organic carbon 
leaching (not measured) 

Equation 1 estimates changes in total below-ground organic 
carbon (TOC) by subtracting respiration and above-ground bio- 
mass production from photosynthetic uptake and includes 
known direct below-ground carbon inputs («), which included 
ploughing in of the original grass root biomass following 
spraying, planted rhizomes and weed roots killed in herbicide 
control. The final term (e) represents carbon exports from the 
site through dissolved organic carbon leaching (DOC) into soil 
water run-off or through soil erosion. Neither of these were 
directly measured at the site and are therefore not included 
although an estimate of the potential level of this export is pre- 
sented in the discussion section below. The carbon content of 
dead plant material is assumed to move either into the soil 
organic carbon (SOC) pool or be respired through decomposi- 
tion as part of ecosystem respiration (Reco). Respiration losses 
are captured in the eddy covariance measurements. Carbon 
contents were assumed to be 44% of Miscanthus biomass and 
42% of grass/weed biomass (J. Clifton-Brown, unpublished). 
This average carbon content for above-ground and below- 
ground Miscanthus biomass combined agrees well with an 
average figure of 46.7% in Beuch et al. (2000) and 45% in Han- 
sen et al. (2004). Equation 2 derives an estimate of the change 
in SOC over time by subtracting the carbon content of below- 
ground biomass production and planted rhizomes. These del- 
tas were compared between years and integrated over the 
three-year study period. Incorporation of carbon from leaf litter 
into the SOC pool was assumed at 26% of the carbon content of 
litter decomposition after 3 years (Hansen et al., 2004), the 
remainder respired back to the atmosphere would be captured 


© 2015 The Authors. Global Change Biology Bioenergy published by John Wiley & Sons Ltd, doi: 10.1111/gcbb.12323 


MISCANTHUS CARBON BALANCE IN LAND-USE CHANGE 7 


in the eddy covariance results. This litter input was assumed to A so 
move from the above-ground biomass pool (ANPP) to the SOC i E E 
pool only at the end of year three in these calculations, 5 5 JU goug RAY 
-= nl = 
although in reality of course this process would have been 2 5 = O ONA tot 
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Results 2 5 
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Si => 
; : ; 5 sg 
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H © oD 
. . ~% [i a fos) o 
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spies g bo = HHHH HHHH 
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Table 3 Interyear comparison of growing season rainfall, soil moisture content (vwc) and mean air/soil temperatures 


Growing season Rainfall Mean vwc Mean air Mean soil 
(May-—October) [mm] 0-15 cm [%] temp. [°C] temp. [°C] 
2012 606.4 46.4 + 3.31 12.6 + 0.83 14.5 + 0.98 
2013 392.2 34.1 + 2.29 13.6 + 1.09 15.5.2 1:28 
2014 311.6 23.3 + 3.91 14.0 + 0.77 15.0 + 0.73 


Soil vwc 0-15 cm (%) 


150 - P) 
€ 
€ 100 ~ 
£ 
£& 504 
© 
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Month 


Fig. 2. (a) Inter-year comparison of monthly average volumetric soil surface water content (vwc, %) and (b) monthly total rainfall 
(mm). Dashed horizontal line in plot (a) shows estimated field capacity, error bars show + SE. 


Baseline biomass 


Following intensive grazing, above-ground biomass was 
minimal across the site so was included in the root core 
sampling and not assessed separately, results suggested 
a baseline biomass stock of 617 + 61.33 g DM m™”. All 
baseline material was killed by spraying with glypho- 
sate in March 2012 and ploughed into the soil in April 
2012. 


Miscanthus scaling parameter 


From the target commercial planting density in 2012 of 
~1.6 rhizomes m7, 1.49 + 0.7 plants m~ were 
counted just prior to harvest in 2013, this had reduced 
to 1.32 + 0.04 plants m~? by the 2014 harvest. There 
were no extensive gaps in the emerging crop with good 
canopy cover, so the second-year reduction suggests 
losses of individual plants distributed randomly across 
the site. These survivorship figures are used to scale 
from plant mean figures to area means for the corre- 
sponding years. 


Biomass production 


2012. The first Miscanthus shoots, after planting on 24th 
April, were recorded on 22nd May with peak canopy 
height for both LUC blocks recorded on 10th October at 
an average across the site of 885.1 + 19.47 mm (+ SE). 
Differences in biomass dry matter production between 
the retained grassland (RG) and the Miscanthus (LUC) 
areas were significant in both above (F = 224.6, P < 0.001) 
and below (F = 4.811, P < 0.05) ground but only signifi- 
cant within each crop type in terms of below-ground RG 
biomass (F = 6.90, P < 0.05) where RGE was lower 
(569.53 + 12.29 g m 7) than RGW (938.91 + 140.13 g m~®. 
Mean above-ground biomass production (weeds and 
Miscanthus) was far less in the LUC area at 
141.49 + 17.15 g m * compared to 834.91 + 45.5 g m7 
in the RG areas. Weeds contributed on average 51% to 
the above-ground LUC biomass and, while this varied 
greatly between quadrats (range between 12 and 81%), 
the differences were not significant between blocks. 
Mean end of year below-ground biomass for the 


LUC blocks was 421.89 + 66.89 gm with 
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754.22 + 85.38 g m~ in the RG blocks. There were no 
significant differences in below-ground biomass, either 
roots or rhizomes, between the LUC blocks. Miscanthus 
rhizomes contributed 121.97 + 17.82 g m” (~25%) to 
the total below-ground LUC biomass with root mate- 
rial (Miscanthus plus weeds combined) contributing 
299.92 + 64.48 g DM m ~. Dividing total bulk rhizome 
weight planted by area indicated a rhizome planting 
weight of 88.07 g DM m giving a rhizome dry mat- 
ter production in the first year of 33.89 + 17.82 g m * 
with root (Miscanthus plus weeds) and rhizome com- 
bined production at 333.82 + 69.23 g m 7. All above- 
ground production in the LUC areas was returned to 
the system, Miscanthus by mowing and leaving on the 
field (1 March 2013) and the weeds killed through her- 
bicide spraying (3rd April). Grass production in the 
RG areas was mown at the same time as the LUC her- 
bicide application and removed from the system. 


2013. Miscanthus growth was far stronger in 2013, first 
shoots noted on 24th April with peak canopy height 
recorded on 3rd October at 2020.19 + 23.35 mm. This 
year there were significant differences in canopy height 
between the two LUC blocks (F = 36.88, P < 0.001) with 
LUCW at 2128.8 + 28.9 mm and LUCE slightly shorter 
at 1911.5 + 21.1 mm. As expected in the maturing crop, 
both above-ground biomass and below-ground biomass 
were far greater in the LUC blocks than in 2012. Total 
above-ground biomass, weeds and Miscanthus together, 
suggested a mean dry matter production of 
793.34 + 68.52 g m ~. The weeds made up far less of the 
total than in 2012, contributing an average of 5.1% to the 
total dry biomass compared to 51% in 2012. Weed contri- 
bution was not significantly different between blocks. 
For Miscanthus alone, peak yield at senescence averaged 
across the site was 731.97 + 64.99 g m”. LUCW pro- 
duced significantly greater (P < 0.01) above-ground dry 
biomass than LUCE with 833.41 + 77.91 g m~? com- 
pared to 630.53 + 87.07 g m™? Results from a second 
round of sampling immediately prior to spring harvest 
in 2014 suggested that of the 750.92 + 67.85 g m™? Mis- 
canthus DM production in 2013 (including in season litter 
drop) 220.2 + 82.39 g m7 were dropped overwinter as 
litter. This represented 29.33% of the total biomass pro- 
duction and left 530.67 + 46.74 g DM m~ or 5.3 Mg 
ha ! remaining as spring harvestable yield from the 
second-year production. 


Below-ground biomass was not significantly different 
for either root or rhizomes between LUC blocks with a 
mean of 607.82 + 60.41 g m”. Mean rhizome weight 
was 459.40 + 54.49 g m`’; subtracting the previous 
year’s mean rhizome weight showed a mean increase in 
dry matter of 337.43 + 57.33 g m * during 2013. Coarse 
roots were measured at 148.42 + 13.24 g m *. 


No significant differences were found between RGE 
and RGW in either above-ground or below-ground 
grassland biomass in 2013 and neither were significant 
differences seen when both blocks were pooled and 
compared between years. The below-ground biomass 
stock was consistent between years with root biomass 
averaging 700.31 + 70.85 g DM m ~ in 2013 compared 
to 754.22 + 85.38 g m ~ in 2012. Above-ground biomass 
proved equally consistent with a mean of 
819.36 + 44.01 g DM m ~ in 2013 compared to the pre- 
vious year (834.91 + 48.27 g m ”). 


2014. For the 2014 growing season sample, quadrats 
were reduced to five and differences between LUC blocks 
were not investigated. Weed production within the eddy 
covariance footprint was minimal in 2014 following effec- 
tive canopy closure and was not measured. Grassland 
production on the RG areas was not monitored. 

First Miscanthus shoots were noted earlier than the pre- 
vious year on Ist April with peak canopy height also 
recorded earlier on 2nd September at 2516.13 + 
39.40 mm. Mean peak yield at senescence across the LUC 
area was 1291.48 + 33.28 g DM m 7. Harvest weights 
after winter were directly measured from the trailer 
weights leaving the field in spring 2015 with moisture 
content assessed through subsampling, off take was 
783.02 g DM m~, suggesting overwinter litter drop and 
harvest losses at 508.46 g DM m`? which represented 
39.4% of the peak biomass production. Below-ground bio- 
mass was measured at 864.44 + 127.63 g DM m ~, of this 
rhizome contributed 716.87 + 102.12 g DM m”, suggest- 
ing an increase of 257.47 + 115.75 g DM m`? during the 
2014 growing season. Coarse root biomass was consistent 
with the previous year at 147.56 + 27.90 g DM m ~. Fig- 
ure 4 shows an interyear comparison of crop canopy 
height development. 


Eddy covariance 


Data collection and retention. Raw data collection at the 
original eddy covariance mast (EC1) was around 90% 
in 2012 and 2013, reducing to 78.2% in 2014 when the 
sensors were returned to the factory during the Mis- 
canthus dormant season for routine calibration. For EC2 
in 2013, this was 92% but reduced to 70.1% during 
2014 due to sensor failure and replacement. Table 4 
details data collection and retention after quality con- 
trol for both masts in all years. Rejection through qual- 
ity control was between 60 and 70% and split evenly 
between automatically rejecting quality control flags of 
2 in the processing output (insufficient turbulence, 
spike detection or non-stationarity), sensor occlusion 
by precipitation and footprint extending beyond the 
ideal fetch. 
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Table 4 Data collection and retention compared between 
masts and years, percentages given relate the maximum possi- 
ble half-hour data points in any year from installation of the 
masts and numbers of data points retained after quality con- 
trol. EC1 and EC2 combined shows the data percentage result- 
ing from a combination of the two EC data sets which was 
then used to drive the gapfilling model 


Total possible Collected % Retained % 


EC1 
2012 17 446 15 126 86.7 5115.0 29.3 
2013 17 520 15:530 88.6 6445.0 36.8 
2014 17 520 13 703 78.2 5652.0 32.3 
EC2 
2013 17 109 15 739 92.0 5973.0 34.9 
2014 17 520 12 280 70.1 5168.0 29.5 
EC1 and EC2 combined Possible Retained % 
2013 17 520 9378 53.5 
2014 17 520 8773 50.1 


Energy balance closure 


Energy balance figures were extremely consistent 
between the two masts; particularly in 2013 (see Fig. 3) 
where fits were 0.65 with R? of 0.9, EC1 in 2012 (not 
shown) was similar with a fit of 0.66 although with 
more scatter in the data with an R? at 0.78. Both 
improved slightly in 2014 with fits of 0.75 and 0.79 for 
EC1 and EC2, respectively. The failure of eddy covari- 
ance to show complete energy balance closure is well 
known in the literature, and this site is no exception, 


and the results here, however, are well within published 
figures. A review of 22 FLUXNET sites revealed a mean 
fit of 0.79 with the lowest at 0.53 with intercepts ranging 
from —32.9 to 36.9 (Wilson et al., 2002). In maize sites, 
Hernandez-Ramirez et al. (2010) calculated a four-year 
mean fit and R? of 0.74 while the grassland to energy 
crop conversions of Zenone ef al. (2013) reported a fit of 
0.70 and R? of 0.80. 


Spatial coverage 


The EC sensors were raised four times during each 
growing season to maintain the target of ~ 2 m above 
mean canopy height (apart from year one, 2012, where 
canopy extension was minimal and sensors were raised 
only once). Figure 4 shows a plot of crop height mea- 
surements with EC sensor height plotted above; from 
this, it can be seen that canopy height was measured 21 
times in 2013 and 16 times in 2014 to monitor the clos- 
ing gap as the canopy developed. In year one the mean 
sensor height above the canopy was 1.91 m (range 
1.6 m to 2.3 m), in year two the mean separation was 
1.94 m (range 1.73 m to 2.25 m), and in year three 
2.48 m (range 1.90 m to 3.1 m). Greater separation in 
subsequent years reflects the increasing roughness of 
the taller crop allowing slightly higher measurements to 
remain within the fetch boundary. Sensors were 
returned to 2 m above stubble height immediately fol- 
lowing harvest in each year. Figure 5 shows the results 
of the footprint model, after quality control, from 2013 
as a georeferenced map of mean bearing and distance 


EC1 2013 EC2 2013 
800 (a) y= 0.65 x + 12.87 (b) y = 0.66 x + 6.67 
600 - R2= 0.882 R2= 0.898 
T -200 ae 
E T T T T T T T ji T 
= -200 0 200 400 600 -200 0 200 400 600 800 
i EC1 2014 EC2 2014 
4 3004 (c) y=0.75 x + 24.16 (a) y=0.79 x + 19.84 
600 + R? = 0.877 R? = 0.896 
400 + 
200 -+ 
0 4 
-200 = T T T T T T T T T 


-200 0 200 400 600 


-200 0 200 400 600 800 


Rn-G (W m~?) 


Fig. 3 Ordinary least square regressions between net radiation minus soil heat storage and sensible plus latent heat. Plots a and b 
compare individual masts (EC1 and EC2) in 2013, c and d compare the two masts in 2014. EC1 in 2012 is not shown but was similar 


to 2013 at y = 0.67x + 23.06 and R? = 0.78. 
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Fig. 4 Line plots showing inter-year comparison of crop canopy height development in each growing season (points and solid line) 
with eddy covariance sensor height adjustment to accommodate the developing crop shown above (dashed line). Sensors were 
returned to 2 m above stubble height immediately following each year’s harvest. 


N 
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Fig. 5 Georeferenced output of the eddy covariance footprint model for 2013 overlain onto the field site, each line represents half 
hour data integration retained after quality control. The direction of the line shows the mean wind vector for each half hour with the 
length indicating the distance that 70% of the estimated flux originated from. Quality control rejects half hour means that extend sig- 
nificantly beyond the ideal fetch area indicated by the measureable fetch boundary marked in the plot. 


for data collection for each half-hour integration. From 
this, it can be seen that coverage from EC1 was largely 
restricted to LUCE, and the addition of EC2 from Jan- 
uary 2013 greatly increased spatial coverage at the site, 
extending the footprint coverage of the measurable fetch 
into LUCW. Where the footprint originated beyond the 
fetch at one mast, data could often be collected at the 


other, providing further quality control constraints were 
met. These data complementarity also extended to peri- 
ods of sensor calibration and other data losses. Consid- 
ering EC1 as the primary data set and gapfilling with 
EC2 improved temporal coverage to 53 and 50% for 
2013 and 2014, respectively (see Table 4). Given that for 
2012 only data from EC1 in LUCE were available, and 
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that EC2 has been used to gapfill EC1 in subsequent 
years, an assessment of sampling uncertainty across the 
site is useful in determining the acceptability of this 
approach. Figure 6 shows the results of the OLS regres- 
sion of paired data points between masts. In 2013, there 
were 2615 paired data points recorded as acceptable 
from both masts in the same half-hour time point, in 
2014 this reduced to 2102. There was close agreement 
between the two towers, fits of 0.98 and 0.99 and inter- 
cepts close to zero show very little systematic error 
within the measurements while the R? scatters suggest a 
sampling uncertainty of 5% in 2013, increasing slightly 
to 7% in 2014. 


Temporal coverage 


Temporal coverage from the single EC1 mast, following 
quality control, was notably poor in 2012. Only 30% of 
the possible maximum half hours were retained to pop- 
ulate the gapfilling model. However, this still repre- 
sented 5115 data points which, as can be seen in Fig. 7a, 
were reasonably well distributed across the year. Data 
coverage improves significantly to just over 50% in 
subsequent years (Fig. 7b,c) with the combination of the 
two data sets. The coastal aspect of the study site 
encouraged good turbulence during night-time hours; 
the site-specific friction velocity (u*) threshold for gap- 
filling and flux partitioning was calculated at 
0.12 m s_!; 27% of total night-time measurements were 
below this compared to 11% of daytime data. Night- 
time u* showed a slightly wider range of values (0.007— 
2.21 ms) compared to daytime (0.009-1.59 m s7’). 
The results of the gapfilling model can be seen in 
Fig. 7(d-f), which clearly shows the stronger develop- 


ment of the crop in 2014 compared to 2013. This strong 
growth tails off earlier in the autumn, however, likely 
the result of the reduced growing season rainfall and 
drying soils (see Fig. 2). 


Net Ecosystem carbon Exchange (NEE) 


Year one. As outlined above, for annual carbon 
exchange and budgets over time years are considered 
nominally from harvest to harvest, in this case March at 
the beginning of one growing season to February in the 
following year. Units are converted from mean flux 
rates (umol CO, m~? s~’) to half-hourly sums (g CO»- 
C m~* hh™’) and integrated into monthly and yearly 
sums. The conversion year saw a net emission of carbon 
from the study site (Fig. 8a); cumulative NEE for the 
year was 200.42 + 9.12 g C m°, which is in close 
agreement with Zenone et al. (2013) who, in three land- 


use conversions from grassland in a similar soil and cli- 
mate, saw an average loss of 241.67 + 27 g C m ~~ dur- 
ing their conversion year. Peak net carbon loss occurred 
during the bare soil cultivation months of April at 
64.94 + 4.61 g C m7 and May at 54.11 + 3.98 g C m”. 
Net uptake of carbon (GPP > Reco) was limited to two 
months: September at —39.98 + 3.29 g C m 7 and Octo- 
ber at —29.08 + 2.42 g C m *. Peak respiration (Reco) 
and uptake (GPP) were both found during the growing 
season of July to October with a maximum monthly 
sum of GPP in September at 178.84 g C m “and Reco in 
August at 174.19 g C m”. 


Year two. Growth was far stronger with a maturing 
crop in the second year, cumulative NEE was 
—61.85 + 14.24 g C m 7. However, despite the net 
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Fig. 6 OLS regressions between paired data points measured at the same time points and retained at both masts after quality control 
and application of the site derived friction velocity threshold. Plot (a) show the paired comparison in 2013 (n = 2615) while plot (b) 
shows 2014 (n = 2102). Dashed line indicates 1 : 1, continuous line indicates model fit. 
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Fig. 7 Fingerprint plots of mean half hourly net ecosystem exchange (NEE) of CO, in pmol CO, m~ s7!, pre (a-c) and post (d-f) 
gapfilling. Plot a shows EC1 data only, b and c show the combined datasets of EC1 and EC2. White areas in plots a—c indicate data 
either missing or rejected in quality control. Colour gradient shows direction and magnitude of CO; exchange, negative values indi- 
cate net uptake to the ecosystem. 
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Fig. 8 Time series plots for years one to three (a-c) of gross photosynthetic uptake [primary productivity (GPP)] and total ecosystem 
respiration (Reco) Shown as monthly sums in grey and black bars with cumulative net exchange of carbon (NEE) shown with the run- 
ning line. For cumulative NEE each plot follows on directly from the one above. NEE error bars indicate propagated standard devia- 
tions from random errors in gapfilling combined with uncertainty estimation from the paired data comparisons. For cumulative NEE 
the x axis indicates the atmosphere/ecosystem boundary, points below this line represent net uptake of carbon from the atmosphere 
(sink) as of that time point and vice versa. GPP and Reco are shown as absolute values. 
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uptake during this year, carrying over NEE from the 
previous year saw the site remaining a net source of 
carbon after the first 2 years with cumulative NEE at 
138.57 + 16.91 g C m ~. Peak net carbon loss was seen 
in December 2013 at 72.84 + 5.24 g Cm”, and the 
site did become a net carbon sink for 4 months 
(August to November 2013) following the strong 
growing season but returned to a net source by the 


end of the year with the accumulation of over-winter 
respiration (Fig. 8b). Maximum GPP and Reo were 
both seen in August at 383.70 and 237.96 g C m, 
respectively. 


Year three. Increased photosynthetic uptake following 
establishment and decreased SOC decomposition saw 
cumulative NEE in year three far exceeding preceding 
years at —755.09 + 35.58 g C m ~. The site returned to 
being a net sink for carbon by the end of June 2014 and 
remained so for the rest of year three (Fig. 8c). Maxi- 
mum net carbon loss was seen during December 2014 at 
66.57 + 5.44 g Cm”. Maximum GPP and Reco were 
seen in July, one month earlier than in 2013, at 501.31 
and 210.67 g C m ”, respectively. Carrying over NEE 
from the previous 2 years saw February 2015 ending 
with cumulative NEE over the 3 years at 
—616.52 + 39.39 g C m ~. 


Carbon budget 


Figure 9 shows the results of Eqns 1 and 2 applied at 
harvest time both to individual years and integrated 
across the entire three-year study period. Year one saw 
an increase in total below-ground carbon (TOC) esti- 
mated at 35.6 g C m°, primarily due to addition of 
ploughed grass material and planted rhizome. However, 
for SOC the mass balance approach suggested a loss of 


143.4 g Cm” over the year. This SOC loss was 
increased during Year 2 with a loss of 374.7 g C m~. 
These losses began to reverse during year 3, TOC gained 
186.8 g C m`? with SOC gaining 73.6 g C m”. Taking 
the 3 years as a whole suggested an increase in TOC, 
including biomass, of 139.26 g C m ~. Litter inputs of 
carbon to the soil were assumed to move from the litter 
pool to the SOC pool at the end of year three, and these 
were calculated at 77.71 g C m ~. This resulted in an 
estimate net loss of SOC after 3 years at 366.8 g C m °. 

Figure 10 shows relative size and change of carbon 
stocks in Mg ha™' across the 3 years; Fig. 10(a) shows a 
close up view of the biomass stock pools over time while 
Fig. 10(b) relates these pools to change in the total carbon 
pool. Total system carbon, including biomass, was gener- 
ally higher than the grazed baseline; particularly at peak 
yield, although error bars (propagated standard errors of 
all individual carbon pools) suggest that this difference 
is unlikely to be significant. Baseline total system carbon 
(above ground and below ground) was estimated at 
81.9 + 3.29 Mg ha! and finished year three (immedi- 
ately prior to harvest) at 83.19 + 3.17 Mg ha '. Results 
from Eqn 2 suggested that SOC specifically had declined 
from a baseline of 78.61 + 3.28 Mg ha! to a low of 
73.43 + 3.07 Mg ha! at the end of year two. This trend 
then reversed during year three with some carbon gains 
and a final SOC estimated at 74.94 + 3.13 Mg ha™!, a 
loss of 4.7% from baseline. 


Discussion 


In this paper, we present a multiyear source/sink 
dynamic of an establishing Miscanthus crop on semi- 
improved grassland and employ a simple mass balance 
approach to estimate changes in soil carbon stocks. 
While the experimental layout, determined by the 
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Fig. 9 Results from Eqns 1 and 2 applied to individual years; units shown are g C m°? (1 g C m~ 


2 


= 0.01 Mg ha7'). Increasing pho- 


tosynthetic uptake (GPP) is seen across the years reflecting the maturing crop. Losses from both total below ground carbon (TOC) 
and soil organic carbon (SOC) are greatest in year two, year three suggests the reversal of the SOC loss trend with small gains. 
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Fig. 10 Carbon stock assessments at peak yield (full senescence) and at spring harvest (following overwinter litter drop). Figures are 
presented in Mg ha", Plot (a) shows a magnified image of individual pools (note y axis scales) while plot (b) shows the size of these 
gains and losses relative to the soil carbon pool and overall carbon stock. Error bars show standard errors propagated across all pools. 


demand for homogeneity of the eddy covariance tech- 
nique, is not statistically ideal, offering only pseudo- 
replication, it does provide measurement at a full 
commercial scale. This is, of course, only a single site 
but the results join a growing body of literature for such 
land-use change with the comparison between two 
eddy covariance masts and the capture of the transition 
period itself offering rare insights. Results given are 
from simple mass balance estimations rather than full 
process models and assumptions are made at several 
points. Unlike rhizomes which are easily identified, sep- 
aration of roots from Miscanthus and grass weeds was 
considered impractical so for year one where grass 
weeds dominated across the site, and all roots found 
were considered to be weeds and assumed to move to 
the SOC pool following herbicide application at the 
beginning of year two. For years two and three, where 
the weed burden was dramatically reduced, all roots 
found were considered to be Miscanthus. There is no 
suggestion that all root material has been captured in 
sampling, fine roots were not recovered so results can 
only be considered a minimum; increases of BNPP in 
Eqn 2 that would follow a full capture of fine root mate- 
rial might suggest slightly higher SOC losses than 
reported. Root turnover rates were not included in the 
calculations, the closely matching root biomass assess- 
ments between years two and three likely reflect the 
inability to capture the fine root production and a figure 
relative to the total root biomass would therefore be 
unreliable. For this mass balance approach, it may not 


be so critical, root stock production in total organic car- 
bon (TOC) is captured through GPP in Eqn 1 but with- 
out a reliable turnover rate cannot be included in Eqn 2. 
Agostini et al. (2015) suggest a root input rate of carbon 
to the SOC pool at between 0.19 and 
0.86 Mg C ha! yr! and discuss further root exudate 
inputs at 0.4 to 1.7 Mg C ha! ae However, they 
point out that measurement data for Miscanthus is lim- 
ited and there may be some suggestion that priming 
effects of root exudates increase respiration through 
mineralization of existing SOC, which could possibly 
outweigh any inputs from root turnover itself. Of 
course, all such respiration would be reflected in the 
Reco figures, suggesting that the TOC results from Eqn 1 
would capture these changes in the total below-ground 
carbon reasonably well, assumptions made in Eqn 2 
about SOC specifically should be treated more as a con- 
servative estimate. 

Litter is taken to be incorporated into the soil at 26% 
of the carbon content of total decomposed litter material 
(i.e. 26% of the carbon content of the difference between 
dropped litter summed across the 3 years and the 
remaining litter found at the beginning of year four). 
This percentage was chosen as a conservative estimate 
taken from the nine-year Miscanthus plantation of Han- 
sen et al. (2004). The alternative, sixteen year plantation, 
in the Hansen et al. (2004) study showed a coefficient of 
retention slightly higher at 29% while Amougou et al. 
(2012) calculated a much higher figure of 35% at their 
five-year-old site. Choosing these slightly higher figures 
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would reduce the estimation of SOC loss after 3 years 
by 2.4 or 7.3%. For the mass balance calculation, litter 
incorporation (i.e. movement from ANPP to SOC) is 
only included at the end of year three. In reality, it 
would have been an ongoing process beginning with 
the decomposition of the first-year crop, which was cut 
and left in the field. Inclusion at this point would have 
the effect of slightly reducing the apparent SOC losses 
in the second year; the respired component of decom- 
posing litter material is captured in the Reco figures 
from the eddy covariance. 

No measurement was made of carbon dissolved in 
soil water run-off (DOC). A comprehensive assessment 
of DOC leaching in Ireland (Kiely et al., 2009) across a 
range of land-uses found losses to be determined by 
pre-existing SOC and precipitation levels. In grasslands 
(n = 32), they found a mean annual export of DOC at 
50.1 + 5.54 kg ha! yr! which agrees well with other 
published studies (Hope et al., 1997; Aitkenhead-Peter- 
son et al., 2007). Grassland SOC levels reported in Kiely 
et al. (2009) ranged from 3.2 to 6.3%, the Aberystwyth 
site falls well within this at 4.8%. Assuming similar 
rainfall, this would increase our estimation of SOC loss 
over the 3 years by 4%, however, rainfall in years two 
and three was particularly low so it might be assumed 
that leaching of carbon in soil water run-off would be 
significantly lower in these years. 

With the discussed caveats, the present work offers 
insights into short-term gains and losses in above- 
ground and below-ground carbon pools with the high- 
est losses of soil carbon seen in year two rather than 
year one. This may be due to a combination of constant 
waterlogging of soils in year one with increasing nutri- 
ent priming of decomposer populations in year two fol- 
lowing weed control and increasing litter drop and root 
exudates (Cheng, 2009; Kuzyakov, 2010). 

The results presented here show that the trend of SOC 
decline was reversed in year three, with the site a consis- 
tent net sink for atmospheric carbon and SOC beginning 
to recoup some of the previous years’ losses. The seques- 
tration figure of 0.7 Mg ha~' for year three is in close 
agreement with the mean C4 carbon accumulation rate 
of 0.78 + 0.19 Mg ha! yr! calculated by Poeplau & 
Don (2014) across 31 European sites. Assuming this 
sequestration rate to be a minimum going forward then 
total recoupment of cultivation losses of SOC would be 
seen by year eight of a potential 15-year crop. This con- 
clusion seems reasonable given results from other stud- 
ies which show increased SOC in grassland conversions 
after 12-16 years (Hansen et al., 2004; Clifton-Brown 
et al., 2007; Schneckenberger & Kuzyakov, 2007) and 
decreases in younger plantations between 3 and 9 years 
old (Schneckenberger & Kuzyakov, 2007; Zimmermann 
et al., 2012). Hansen et al. (2004) and Zatta et al. (2013) 


showed no change or a small insignificant decrease after 
9 and 6 years, respectively. Finally, to summarize in 
terms of potential fossil fuel substitution, combined har- 
vest offtake from years two and three was 
13.1 Mg DM ha™. With energy intensity of Miscanthus 
estimated between 15 and 17.5 MJ (kg DM)! (Hastings 
et al., 2008; Felten et al., 2013) and coal at 24 MJ kg! 
each hectare could have potentially replaced 8.2-9.6 Mg 
of coal at a co-fired power station by the end of year 
three. Assuming a typical carbon content of coal at 80%, 
this would offset 6.6-7.6 Mg of carbon emission, more 
than double the estimated loss of SOC over the first three 
years of establishment. 
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